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The LEP experiments have measured event shapes using data taken at e + e — center-of-mass energies 
y/s ranging from 91 GeV to 209 GeV. Using the final LEP event shape measurements, a combined 
value of the strong coupling constant a s (Mz) has been extracted. Events with photon radiation have 
been used to extend the measurements to lower center-of-mass energies \fs' to study the running of 
a s . Alternative a B measurements using four-jet rates have also been performed. 



1 Data Sample 

Between 1989 and 2000 the four LEP exper- 
iments have taken data at e + e~ center-of- 
mass energies y/s ranging from 91 GeV to 209 
GeV. Each experiment has collected several 
million multi-hadronic events at y/s = Mz- 
At higher energies about 10 3 events have 
been recorded per energy point. After a basic 
event selection, no significant background re- 
mains at y/s = Mz- In the range y/s > 2A/\y 
the background rate is about 10 — 15% due to 
four-fermion production in e + e _ — > W + W~ 
events. For an accurate determination of y/s 
it is also important to identify and recon- 
struct events with initial and final state pho- 
ton radiation. 

2 Event Shape Observables 

The properties of hadronic events may be 
described by a set of event shape observ- 
ables. These may be used to characterize 
the distribution of particles in an event as 
"pencil-like", planar, spherical, etc. They 
can be computed either using the measured 
charged particles and calorimeter clusters, or 
using the hadrons or partons in simulated 
events. In order that predictions of pertur- 
bative QCD be reliable, it is necessary that 
the value of the observable be infra-red stable 
(i.e. unaltered under the emission of soft glu- 
ons) and collinear stable (i.e. unaltered under 
collinear parton branchings). The following 
event shapes are considered here: 



Thrust T: defined by the expression 

« V EJm 

where pi is the three-momentum of par- 
ticle i. The thrust axis Ht is the direc- 
tion ft which maximizes the expression 
in parentheses. A plane through the ori- 
gin and perpendicular to ftr divides the 
event into two hemispheres. 

C-parameter: The linearized momentum 
tensor is given by 

e a/3 = J2MPi )/\Pi\ (o ^ = lj2 ,3), 
Ei ml 

where the sum runs over particles, i, and 
a, [3 denote the Cartesian coordinates of 
the momentum vector. The three eigen- 
values Xj of this tensor define C through 

C — 3(AiA2 + A2A3 + A3A1). 

Heavy Jet Mass M h- The hemisphere in- 
variant masses are calculated using the 
particles in the two hemispheres. Mh is 
the heavier of the two masses. 

Jet Broadening and B\y'- These are 
defined by computing the quantity 

'Eieij, \Pi x ft T \ 



B, 
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for each of the two hemispheres. The two 
observables arc defined by 

Bt = B1+B2 and Bw = max(.Bi, B2), 
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Figure 1. (left) 1 - T measured by OPAL in the range = 91 - 197 GeV. (center) M H measured by DELPHI 
at y/s = 207 GeV. The subtracted four-fermion background is shown separately. The OPAL and DELPHI 
data are compared to PYTHIA (JETSET), HERWIG and ARIADNE MC simulations, (right) B T measured 
by ALEPH in the range ^fs = 91 — 206 GeV. The continuous line shows the region where the C(c«g)+NLLA 
fit has been performed. 



where Bt is the total and B\y is the wide 
jet broadening. 

Transition between 2- and 3-jets (1/23) : 

The value of the jet resolution param- 
eter, y cu t, at which the event makes a 
transition between a 2-jet and a 3-jet 
assignment, for the Durham jet finding 
scheme . 

The final results on event shapes have re- 
cently been published by all four LEP col- 
laborations 2 - 3 ' 4 ' 5 . Examples are shown in 
Fig. 1. The agreement between data and 
Monte Carlo (MC) simulations is in general 
very good which justifies the use of simple 
bin- by-bin efficiency corrections. 

The QCD calculations for the cumula- 
tive cross-section R(y) = ^-^dy as a 
function of the event shape variable y{y — 
1 — T, C, Mh, Bt, B\y, j/23) can be performed 
in the next-to-leading-logarithmic approxi- 
mation (NLLA) or using a fixed order 0(a?) 
calculation. The 0(a%) prediction is ex- 
pected to be best for the multi-jet region 
(high y), whereas the NLLA prediction is 
best for the two-jet region (low y). The 
"log(i?)" matching scheme is adopted for 
combining the 0{a%) and NLLA predictions. 



In the log(i?) matching scheme the terms up 
to O{ol) in the NLLA expression are re- 
placed by the 0(a%) terms from logi?e>( Q 2). 
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Figure 2. Combined a s (Mz) for the six event shape 
variables measured by the four LEP experiments. 
The inner error bar is the statistical, the outer er- 
ror bar is the total uncertainty 6 . 



Each experiment has performed a s fits 
using these predictions for the six event shape 
variables at different i/s. A preliminary com- 
bination of these results, taking into account 
the correlations between the sources of un- 
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certainties, yields 6 
a B (M z ) = 
0.1202 ± 0.0003(stat) ± 0.0007(expt) 
±0.0015(hadr) ± 0.0044(theo) 

with x 2 /d.o.f. = 93/166. The combined a s 
measurements for each event shape variable 
are shown in Fig. 2 and the combined a s mea- 
surements for each experiment in Fig. 3. 
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Figure 3. Combined a B (Mx) for the four LEP exper- 
iments using the measurements from the six event 
shape variables. The inner error bar is the statisti- 
cal, the outer error bar is the total uncertainty 6 . 
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Figure 4. The running of a B as a function of Q = \/s 
using the combined LEP data. The inner error bar is 
the statistical uncertainty, the outer error bar is the 
total uncertainty 6 . 



ADNE, PYTHIA. The theoretical uncer- 
tainty includes the variation of the renormal- 
ization scale, 0.5 < x^ < 2, the variation of 
the logarithmic rescaling factor, 2/3 < xl < 
3/2, and of kinematic cut-off parameters. In 
addition, the log(i?) matching scheme is re- 
placed by the R matching scheme 7 . 

The running of a s using the combined 
data is shown in Fig. 4. Although the statis- 
tical precision of the LEP1 data at ^fs = Mz) 
is much higher, the LEP1 and LEP2 data 
have about equal weight in the combined fit 
because of the smaller hadronization and the- 
ory uncertainties at larger ^/s. 

3 a s from Radiative Events 

Radiative hadronic events have been used to 
extend the range of center-of-mass energies 
below yfs = Mz- Assuming that initial state 
or final state photons do not interfere with 
the QCD process 8 , a measurement of a s at a 
reduced center-of-mass energy y / s' is possible 
by identifying isolated high energy photons in 
the detector 9 . 
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Figure 5. The running of a s as a function of Q = 
VV using radiative hadronic events measured by 
OPAL 10 . The result of the a s fit is shown as contin- 
uous line. 



The experimental uncertainty includes 
uncertainties from cut variation, detector 
corrections and background subtraction. The 
hadronization uncertainty is determined by 
performing the hadronization correction us- 
ing different MC generators (HERWIG, ARI- 



In the analysis performed by OPAL 10 
photons have been separated from the dom- 
inant 7T° background by using a likelihood 
based on cluster shape fits in the calorime- 
ter. The running of a s as function of v / s' is 
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Table 1. Summary of the a B (Mz) results presented in this talk. 





data 


theory 


a s (Mz) stat expt hadr theo 


LEP 
OPAL 
OPAL 
DELPHI 


event shapes 
radiative events 

i?4 /Durham 
i?4/ Cambridge 


£>(a 2 )+NLLA 
0(a s 2 )+NLLA 
0(c^)+NLLA 
0( a 3) +x opt 


0.1202 0.0003 0.0007 0.0015 0.0044 
0.1176 0.0012 ±°;°°fg 
0.1208 0.0006 0.0021 0.0019 0.0024 
0.1175 0.0005 0.0010 0.0027 0.0007 



shown in Fig. 5. The result of the a s fit is 
given in Table 1. 

4 Power Law Corrections 

Non-perturbative effects in event shape ob- 
servables are usually suppressed by powers 
of 1/Q 11 . The corresponding hadronization 
corrections at LEP energies are of the order 
10%. These hadronization corrections can be 
determined using MC generators or by ana- 
lytical power law calculations. These calcu- 
lations introduce one additional phenomeno- 
logical parameter ao, 

1 



Mi Jo 



a s (k)dk, 



which measures the effective strength of the 
strong coupling up to an infrared matching 
scale of jii w 1 GeV. The parameter ao is 
expected to be universal and must be deter- 
mined by experiment. 
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Figure 6. Contours of confidence levels for simulta- 
neous measurements of a B and ao measured by L3 
using the first moments of event shape variable. The 
unweighted averages are shown as shaded bands. A 
4-jet observable (D parameter) is also shown. 
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Figure 7. Contours of confidence levels for simultane- 
ous measurements of a B and ao compared to the a B 
measurement using MC corrections (shaded band). 



The result of such a combined fit, per- 
formed by L3 4 for the first moments of the 
event shape variables, is shown in Fig. 6. 
The six values of ao obtained from the event- 
shape variables do not agree well. A similar 
result has been obtained by DELPHI 3 . A 
large spread of a values has also been ob- 
served by ALEPH 2 using event shape dis- 
tributions. The resulting a s is significantly 
lower than the result obtained using MC 
hadronization corrections (Fig. 7). 

5 a s from Four-jet Rates 

An alternative method to determine a s is 
to measure four-jet rates, since the jet rates 
are predicted as functions of the jet resolu- 
tion parameter with a s as free parameter. 
OPAL 12 has measured the four-jet rates us- 
ing the Durham jet finding algorithm 1 as 
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a function of the jet resolution parameter 
(Fig. 8). The strong coupling a s is fitted us- 
ing a 0(a!s)+NLLA matched calculation of 
the four-jet rates. 
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Figure 8. A fit to the four-jet rate Ra measured at 
y/s = 198 GeV. 
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Figure 9. Fits to the four-jet rates R4 measured at 
y/s = Mz using the Cambridge algorithm. The grey 
band shows a fit with the physical scale (x M = 1) and 
the cross-hatched band a fit with optimized scales. 



In a similar analysis, DELPHI 13 has 
measured the four-jet rates using the Durham 
and the Cambridge 14 algorithms. The data 
are fitted using a fixed order O(a^) calcu- 
lation with both a s and the renormalization 
scale free parameter. To determine the 

experimentally optimized scale, a two param- 
eter fit with a s and x° pt as free parameters 
has also been performed (Fig. 9). Due to a 
smaller theoretical uncertainty the a s fit us- 
ing the Cambridge algorithm is quoted as fi- 
nal result (Table 1). 

6 Summary 

A preliminary a s combination using the final 
event shape measurements of all four LEP 
experiments has been presented. The uni- 
versality of power law corrections has been 
studied. Radiative events and four-jet events 
have been used to study the running of a s . 
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